Thermalization of a nonequilibrium electron-positron-photon plasma 
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Starting from a nonequilibrium configuration we analyse the essential role of the direct and the 
inverse binary and triple interactions in reaching an asymptotic thermal equilibrium in a homoge- 
neous isotropic electron-positron-photon plasma. We focus on energies in the range 0.1-10 MeV. 
We numerically integrate the integro-partial differential relativistic Boltzmann equation with the 
exact QED coUisional integrals taking into account all binary and triple interactions in the plasma. 
We show that first, when detailed balance is reached for all binary interactions on a timescale 
tk ^ 10"^'* sec, photons and electron-positron pairs establish kinetic equilibrium. Successively, when 
triple interactions fulfill the detailed balance on a timescale teq < 10~^^sec, the plasma reaches 
thermal equilibrium. It is shown that neglecting the inverse triple interactions prevents reaching 
thermal equilibrium. Our results obtained in the theoretical physics domain also find application in 
astrophysics and cosmology. 

PACS numbers: 52.27.Ep; 05.20.Dd 
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An electron-positron plasma is of interest in many 
fields of physics and astrophysics: the early universe p|, 
gamma-ray bursts B, active galactic nuclei 0], the cen- 
ter of our Galaxy hypothetical quark stars [1] and 
ultraintense lasers [6|. 

A detailed study of the relevant processes and possi- 
ble equilibrium configurations in an optically thin pair 
plasma are given in [7[. In all the above-mentioned ap- 
plications the precise knowledge of the optically thick 
plasma evolution is required. In this case there exists 
only a qualitative description and an assumption of ther- 
mal equilibrium is often adopted without explicit proof 

In this Letter we consider a uniform isotropic electron- 
positron-photon plasma in the absence of external elec- 
tromagnetic fields and we describe its evolution start- 
ing from arbitrary nonequilibrium initial conditions up 
to reaching thermal equilibrium. We are interested in 
the range of final temperatures in thermal equilibrium, 
bracketing the electron rest mass energy 



0.1MeV<rt,, < 10 MeV. 



(1) 



These boundaries are required for the study of electron- 
positron pairs in absence of the production of other par- 
ticles such as muons. We assume that the energy den- 
sity of the plasma is constant and is, correspondingly, 
in the range 1.6 lO^^-^ < p < 3.8 10^°-^. The rel- 
ative number densities at thermal equilibrium will be 
3.1 lO^^cm-^ < nth < 7.9 lO^^cm^^. 

We adopt a kinetic description for the distribution 
functions of electrons, positrons and photons. In our case 
the plasma parameter is small, g — (nr^)^^ 1, where 
r£) is the Debye length, and therefore we use one-particle 
distribution functions. Besides, in our case electrons and 
positrons are non-degenerate. We solve numerically the 
relativistic Boltzmann equations Q which for homoge- 
neous and isotropic distribution functions of electrons, 



positrons and photons reduce to 
c dt ^ 

9 



(2) 



where fi{e,t) are their distribution functions, the index i 
denotes the type of particle, e is their energy, and 77* and 
Xi are the emission and the absorption coefficients for 
the production of a particle of type "i" via the physical 
process labeled by q. 

In order to solve equations ^ we use a finite difference 
method by introducing a computational grid in the phase 
space to represent the distribution functions and to com- 
pute coUisional integrals ^ . The result of this procedure 
is the stiff system of ordinary differential equations to be 
solved with the implicit Gear method [lo[. For binary 
interactions we use exact QED matrix elements [ll]| . For 
triple interactions we compute emission and absorption 
coefficients following Svensson [12[. The Compton scat- 
tering of photons, for instance, is described by 0] 



where 

U'k'p';k,p 



J dk'dpdp'wk'p';k,p/7(k', <)/±(p', t), 
y"(ik'dpdp'-u;k'p';k,p/7(k, t)f±{p, t), 



(3) 
(4) 

(5) 



is the corresponding transition probability, k = (e^/c, k) 
and p — (fe/c, p) are four-momenta of photon and 
positron (electron), primes denote particles after the in- 
teraction, and Mfi is the matrix element for the consid- 
ered process. 

For such a dense plasma coUisional integrals in ^ 
should include not only binary interactions, having order 
a in Feynmann diagrams, where a is the fine structure 



2 



Binary interactions 


Radiative variants 


M0ller, Bhabha 

± ±1 ±" ±"' 
e e e e 


Bremsstrahlung 

e±e±'^e±"e±"'7 
e^e^^e^'e^'7 


Single Compton 


Double Compton 

e^7<-^e^'7'7" 


Pair production 
and annihilation 

77'<— >e^e^ 


Radiative pair production 
and three photon annihilation 

77'^e^e^7" 
e 7^e e^e 
e^e^<— >77'7" 



TABLE I; Microphysical processes in the pair plasma. 



constant, but also triple ones, having order ll|. We 
consider all possible binary and triple interactions be- 
tween electrons, positrons and photons as summarized in 
Tab.H 

Each of the above reactions is characterized by the 
corresponding timescale and optical depth. For Compton 
scattering, for instance, we have 



tc 



1 



aTn±c 



aTn±Ro, 



(6) 



where ax is the Thomson cross-section and n± is the 
number density of pairs. There are two timescales in 
our problem that characterize the condition of detailed 
balance between direct and inverse reactions, ^ tcs for 
binary and a~^tca for triple interactions respectively. In 
the first phase of the system evolution the binary inter- 
actions are found to have a predominant role. Starting 
from arbitrary distribution functions we find a common 
development: at the time tcs the distribution functions 
always have evolved in a functional form on the entire 
energy range, depending only on two parameters. We 
find in fact for the distribution functions the expressions 



f^{e) = exp 



(7) 



with chemical potential ipi = ^'^^i and temperature 
6i = where e = is the energy of the particle, 

TOg is the electron mass and /cb is Boltzmann's constant. 
Such a configuration corresponds to a kinetic equilibrium 
[ist , 14 1 in which electrons, positrons and photons ac- 
quire a common temperature and nonzero chemical po- 
tential. At the same time we found that triple interac- 
tions become essential for t > tcs, after the establishment 
of kinetic equilibrium. Such triple interactions, both di- 
rect and inverse, are indeed essential in achieving the 
thermal equilibrium. 

In ^ analogously to the temperature, defining the av- 
erage kinetic energy in the system, the chemical poten- 
tial represents deviation from the thermal equilibrium 




10"'' 10" 10"'' 10"'' 10"" 10"" 10"' 

t,s 



FIG. 1: Depencence on time of concentrations of pairs (black), 
photons (red) and both (thick) when all interactions take 
place (solid). Upper (lower) figure corresponds to the case 
when initially there are mainly pairs (photons) . Dotted curves 
on the upper figure show concentrations when inverse triple 
interactions are neglected. In this case an enhancement of 
the pairs occurs with the corresponding increase in photon 
number and thermal equilibrium is never reached. 



through the relation ip = 9hi{n/neq), where n^q are con- 
centrations of particles in thermal equilibrium. We do 
not absorb the chemical potentials into the normaliza- 
tion factors since they depend on time and describe the 
approach to thermal equilibrium. 

The results of numerical simulations are reported be- 
low. We choose two limiting initial conditions with flat 
spectra: a) electron-positron pairs with a 10^^ energy 
fraction of photons and b) the reverse case, i.e., pho- 
tons with a 10~^ energy fraction of pairs. Our grid 
consists of 60 energy intervals and 16 x 32 intervals for 
two angles characterizing the direction of the particle 
momenta. In both cases the total energy density is 
p = 10'^'*erg/cm^. In the first case initial concentration 
of pairs is 3.1 lO^^cm^'^, in the second case the concen- 
tration of photons is 7.2 lO^^cm"^. 

In Fig. [T] we show concentrations of photons and pairs 
as well as their sum for both our initial conditions. After 
calculations begin, concentrations and energy density of 
photons (pairs) increase rapidly with time, due to anni- 
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Interaction 


Parameters of distribution functions 


Compton scattering 




Pair production 


Lp^ — <^-|-, if 0^ — 9± 


Tripe interactions 


if^, ip± ^ 0, if 6^1 = 0± 



TABLE 11: Relations between parameters of equilibrium DFs 
fulfilling detailed balance conditions for each of the reactions 
shown in Tab. |I] 



hilation (creation) of pairs by the reaction 77'<->e^e^. 
Then, in the kinetic equilibrium phase, concentrations of 
each component stay almost constant, and the sum of 
concentrations of photons and pairs remains unchanged. 
Finally, both individual components and their sum reach 
stationary values. If one compares and contrasts both 
cases as reproduced in Fig. [T] one can see that, although 
the initial conditions are drastically different, in both 
cases the same asymptotic values of the concentration 
are reached. 

We now describe in detail the case when initially pairs 
dominate. One can see in Fig. [2] that the spectral density 
of photons and pairs 



(8) 



1 for 

15„ 



where (3± = \J\ ~ (mgC^/e)^ for pairs and fi-y — 
photons, can be fitted already at tk ~ 20ics — 7 10~^''sec 
by distribution functions ([7]) with definite values of tem- 
perature Ok{tcs) ~ 1-2 and chemical potential ipkitk) ~ 
—4.5, common for pairs and photons. As expected, after 
tk the distribution functions preserve their form ([7]) with 
the values of temperature and chemical potential chang- 
ing in time, as shown in Fig. [31 As one can see from Fig. [3] 
the chemical potential evolves with time and reaches zero 
at the moment tth ~ ct^^tk — 7 lO^^'^sec, corresponding 
to the final stationary solution. 

We now discuss the results. Let us consider the distri- 
bution functions ([7]) with different temperatures 9i and 
chemical potentials ipi for pairs and photons. The re- 
quirement of vanishing reaction rate for the Compton 
scattering f±f^ = f^fj leads to the equal temperature 



of pairs and photons 9± = = 9k, see also [13[,[1J|. In 
this way the detailed balance between any direct and the 
corresponding inverse reactions shown in Tab. U leads to 
relations between 9 and (p collected in Tab. |TT1 

These relations are not imposed, but are verified 
through the numerical calculations. This is a powerful 
tool to verify the consistency of our approach and nu- 
merical calculations. These relations were obtained for 



the first time, to our knowledge, in 1J| and then later in 
[ist for binary reactions. 

From Tab. lU one can see that the necessary condition 
for thermal equilibrium in the pair plasma is detailed 
balance between direct and inverse triple interactions. 
This point is usually neglected in the literature where 




e -m c , keV 

r e 



FIG. 2: Spectra of pairs (upper figure) and photons (lower 
figure) when initially only pairs are present. The black curve 
represents the results of numerical calculations obtained suc- 
cessively a.t t = 0, t = tk and t = tth (see the text). Both 
spectra of photons and pairs are initially taken to be flat. 
The yellow curves indicate the spectra obtained form ((7} at 
t — tk- The perfect fit of the two curves is most evident in the 
entire energy range leading to the first determination of the 
temperature and chemical potential both for pairs and pho- 
tons. The orange curves indicate the final spectra as thermal 
equilibrium is reached. 



there are claims that the thermal equilibrium may be es- 
tablished with only binary interactions 15]. In order to 
demonstrate it explicitly we also show in Fig. [1] the de- 
pendence of concentrations of pairs and photons when 
inverse triple interactions are artificially switched off. In 
this case, see dotted curves in the upper Fig. [U after 
kinetic equilibrium is reached concentrations of pairs de- 
crease monotonically with time, and thermal equilibrium 
is never reached. 

The existence of a non-null chemical potential for pho- 
tons indicates the departure of the distribution function 
from the one corresponding to thermal equilibrium. Neg- 
ative (positive) value of the chemical potential generates 
an increase (decrease) of the number of particles in or- 
der to approach the one corresponding to the thermal 
equilibrium state. Then, since the total number of par- 
ticles increases (decreases), the energy is shared between 
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positrons and photons, the sufficient condition to obtain 
an isothermal distribution on a causally disconnected 
spatial scale R > ctth — lO^^cm is the request of con- 
stancy of the energy density on such a scale as well as, 
of course, the invariance of the physical laws. 

We have considered the evolution of an initially 
nonequilibrium optically thick electron-positron-photon 
plasma up to reaching thermal equilibrium. Starting 
from arbitrary initial conditions we obtain kinetic equi- 
librium from first principles, directly solving the relativis- 
tic Boltzmann equation with coUisional integrals com- 
puted from QED matrix elements. We have demon- 
strated the essential role of direct and inverse triple in- 
teractions in reaching thermal equilibrium. Our results 
can be applied in the theories of the early universe and 
of gamma-ray bursts, where thermal equilibrium is pos- 
tulated at the very early stages. These results can in 
principle be tested in laboratory experiments in the gen- 
eration of electron-positron pairs. 

We thank the anonymous referee for comments which 
have improved the comprehension of our results. 



FIG. 3: Time dependence of temperatures, measured on the 
left axis (solid), and chemical potentials, measured on the 
right axis (dotted), of electrons (black) and photons (red). 
The dashed lines correspond to the reaching of the kinetic 
(~ 10~^'*sec) and the thermal (~ 10~^^sec) equilibria. Upper 
(lower) figure corresponds to the case when initially there are 
mainly pairs (photons). 



more (less) particles and the temperature decreases (in- 
creases), see fig. 131 Clearly, as in thermal equilibrium is 
approached, the chemical potential of photons is zero. 

In our example with the energy density 10^'*erg/cm^ 
the thermal equilibrium is reached at ~ 7 10~^'^sec with 
the final temperature Tth ~ 0.26 MeV. For a larger 
energy density the duration of the kinetic equilibrium 
phase, as well as of the thermalization timescale, is 
smaller. In our entire temperature range ([T]) we deal 
with a non-degenerate plasma. 

Our results, obtained for the case of an uniform 
plasma, can only be adopted for a description of a phys- 
ical system with dimensions Rq ^ = 4.3 lO^^cm. 

The assumption of the constancy of the energy density 
is only valid if the dynamical timescale tj^y^ — {^^) 
of the plasma is much larger than the above timescale 
tth which is indeed true in all the cases of astrophysical 
interest. 

Since we get thermal equilibrium already on the 
timescale tth ^ 10~^^sec, and such a state is indepen- 
dent of the initial distribution functions for electrons, 
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